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Abstract. This article consists of two parts. In the first part, we review the most recent 
proofs estabhshing quadratic Morawetz inequalities for the nonlinear Schrodinger equation 
(NLS). We also describe the applications of these estimates to the problem of quantum 
scattering. In the second part, we generalize some of the methods developed for the NLS 
by many authors to the case of Gross-Pitaevskii (GP) hierarchies. In particular, we prove 
both regular and interaction Morawetz identities for the GP hierarchy, which appear here 
for the first time in the literature. 

1. Introduction 

In this paper, we discuss some recent a priori estimates for the solutions to the power 
type nonhnear Schrodinger equation 

' ' 1.1 

u(x,0) = uo{x) € H'{W) 

for any p > 1 and A € M. This involves the proof of Morawetz identities and Morawetz 
inequalities for the NLS. Moreover, we derive generalizations of both regular (one-particle), 
and interaction Morawetz identities for Gross-Pitaevskii hierarchies. 

Gross-Pitaevskii (GP) hierarchies are infinite systems of coupled linear non-homogeneous 
PDEs that appear naturally in the derivation of the nonlinear Schrodinger equation, as the 
dynamical mean field limit of the manybody quantum dynamics of Bose gases (under the 
Gross-Pitaevskii scaling). Let q G {2,4} and A G {—1,1}. The q-GP hierarchy on is 
an infinite system of coupled PDEs for functions (referred to as marginal density matrices) 
7(^)(t,xi, . . . ,Xk,x[, . . . ,x'f^) for k eN, t and Xi,x[ G R'^ for i e {1, . . . ,k}: 

ida^'^ = -E (^-. - ^-;) ^^'^ + ^^^+§7^'+'^ (1-2) 
i=i 

where the inhomogeneous term is given by 

^^..P"**' = E (B^.« - B-,^, ,,,.<^-i)) . (1,3) 
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accounting for the interactions between 1 + | particles. The marginal density matrices are 



hermitean, 7^'^H3ifc'i£fc) = T'-'^-'fefci^^fc)) ^-^^ satisfy the property of admissibility, that is, 

-f^^\t,xi,...,Xk;x[,...,x'^) = j dy-f^'''^^\t,xi,...,Xk,y;x[,...,x'k,y) (1.4) 
for all fc G N, with normalization J dx^^^\t,x]x) = 1. 

1.1. Nonlinear Schrodinger equation. Equation (II. ip is the Euler-Lagrange equation 
for the Lagrangian density 

L(n) = -^A(H2) + A^|nri. 

Several invariants of the equation lead to conservation laws that play a central role in the 
well-posedness theory of the problem. The homogeneous Sobolev norm H'^" is invariant 

2_ 2 

under the scaling u^{x,t) = fi j^), when Sc = § — This critical regularity 

already partitions the general problem into different sub-problems with a varied degree of 
difficulty. Two important cases that are well understood correspond to the energy subcrit- 
ical problem for 1 < p < 1 + n > 3, 1 < p < oo, n = 1, 2 and initial data in H^, and 
the subcritical problem forl<p<l-|-^, n > I and initial data. Time translation 
invariance leads to energy conservation 

E(u)(t) = - [ \Vu(t)\^dx + [ \u(t)\P+^dx = E(uo). (1.5) 

2 J P+l J 

In the defocusing case, (A = 1), this law provides an a priori bound that extends any 

local energy subcritical solution to an interval of arbitrary length. The equation is also 

invariant under phase rotation which leads to the conservation of mass law 

\W{t)\\L2 = iimo||l2. (1.6) 

This implies for any A the existence of global-in-time solutions for solutions evolving from 
finite mass data in the L^-subcritical case of p < 1 + -. For the L^— critical problem, 
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p = 1 + ^, the time of the local theory depends on the profile of the initial data, in addition 
to their norm, and thus the conservation law is not immediately useful. Space translation 
invariance leads to momentum conservation 

p{t) = 3 / uVudx, (1.7) 

a quantity that has no definite sign. It turns out that one can also use this conservation 
law in the defocusing case and prove monotonicity formulas that are very useful in studying 
the global-in-time properties of the solutions at t = oo. For most of these classical results 
the reader can consult [3], [13] . 

The study of the problem at infinity is an attempt to describe and classify the asymp- 
totic behavior-in-time for the global solutions. To handle this issue, one tries to compare 
the given nonlinear dynamics with suitably chosen simpler asymptotic dynamics. For the 
semilinear problem (jl.ip . the first obvious candidate for the simplified asymptotic behavior 
is the free dynamics generated by the group S{t) = e~**^. The comparison between the two 
dynamics gives rise to the questions of the existence of wave operators and of the asymp- 
totic completeness of the solutions. In the literature nowadays, the problem of asymptotic 
completeness is synonymous with the problem of scattering. More precisely, we have: 

i) Let v^{t) = S{t)u^ be the solution of the free equation. Does there exist a solution u 
of equation (jl.ip which behaves asymptotically as u+ as t oo, typically in the sense that 
\\u{t) — v+Wf^i 0, as t — 7- oo. If this is true, then one can define the map 0+ : n(0). 
The map is called the wave operator and the problem of existence of u for given u-^ is 
referred to as the problem of the existence of the wave operator. The analogous problem 
arises as t ^ — oo. 

ii) Conversely, given a solution u of (jl.ip . does there exist an asymptotic state n+ such 
that = S{t)u+ behaves asymptotically as u{t), in the above sense. If that is the case 
for any u with initial data in X for some n+ E X, one says that asymptotic completeness 
holds in X. 

In effect the existence of wave operators asks for the construction of global solutions that 
behave asymptotically as the solution of the free Schrodinger equation while the asymptotic 
completeness requires all solutions to behave asymptotically in this manner. It is thus not 
coincidental that asymptotic completeness is a much harder problem than the existence of 
the wave operators (except in the case of small data theory which follows from the iterative 
methods of the local well-posedness theory). 

Asymptotic completeness for large data not only require a repulsive nonlinearity but 
also some decay for the nonlinear solutions. In the ideal scenario the nonlinear solution will 
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obey the same decay properties of the linear Schrodinger equation. This decay of the hnear 
problem can already establish weak quantum scattering in the energy space for example, 
but to say more, we usually need decay in an norm for the nonlinear solution. In this 
paper, we provide a summary of recent results that demonstrates a straightforward method 
to obtain such estimates by taking advantage of the momentum conservation law (|1.7p . A 
key example of these ideas is contained in the following generalized virial inequality of Lin 
and Strauss [33] • (We recall the proof of this inequality in the next section.) 



/ {-AAa{x))\u{x,t)\^dxdt + X I {Aa{x))\u{x,t)\P+^dxdt < C (1.8) 

where a{x) is a convex function, tt is a solution to (jl.ip . and C a constant that depends 
only on the energy and mass bounds. 

An inequality of this form, which we will call a one-particle inequality, was first derived 
in the context of the Klein-Gordon equation by Morawetz |34j . and then extended to the 
NLS equation in [33j. The inequality was applied to prove asymptotic completeness first 
for the nonlinear Klein-Gordon and then for the NLS equation in the papers by Morawetz 
and Strauss, p5j, and by Lin and Strauss, [33] for slightly more regular solutions in space 
dimension n > 3. The case of general finite energy solutions for n > 3 was treated in j25j 
for the NLS and in [23] for the Hartree equation. The treatment was then improved to the 
more difficult case of low dimensions by Nakanishi, |36] . [37] . 

The bilinear a priori estimates that we outline here give stronger bounds on the solutions 
and in addition simplify the proofs of the results in the papers cited above. For a detailed 
summary of the method see [24J. In the original paper by Morawetz, the weight function 
that was used was a{x) = \x\. This choice has the advantage that the distribution — AA(|^) 
is positive for n > 3. More precisely it is easy to compute that Aa{x) = and that 

{87r(5(x), if n = 3 
^^^^^f^, if ri>4. 

In particular, the computation in (jl.Sp gives the following estimate for n = 3 and A positive 

\u{t,0)\^dt+ \ I dxdt<C. (1.9) 

Jr Jr^xr fI 

Similar estimates are true in higher dimensions. The second, nonlinear term, or certain local 
versions of it, have played central role in the scattering theory for the nonlinear Schrodinger 
equation, [2], [25], [26], [33]- The fact that in 3d, the bi-harmonic operator acting on the 
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weight a(x) produces the 5— measure can be exploited further. In [13], a quadratic Morawetz 
inequahty was proved by correlating two nonlinear densities pi{x) and P2{y) and define as 
a{x, y) the distance between x and y in 3d. The authors obtained an a priori estimate of 
the form /j^Sxr \u{x,t)\'^dx < C for solutions that stay in the energy space. A frequency 
localized version of this estimate has been successfully implemented to remove the radial 
assumption of Bourgain, [2], and prove global well-posedness and scattering for the energy- 
critical (quintic) equation in 3d, [Mj. For n > 4 new quadratic Morawetz estimates were 
given in 05]. Finally in |10| and in [38] these estimates were extended to all dimensions. 

We should mention that taking as the weight function the distance between two points 
in M" is not the only approach, see [11] for a recent example. Nowadays it is well under- 
stood that the bilinear Morawetz inequalities provide a unified approach for proving energy 
scattering for sub-critical solutions of the NLS when p > 1 + ^ {L? super-critical nonlin- 
earities). This last statement has been rigorously formalized only recently due to the work 
of the aforementioned authors, and a general exposition has been published in |24j . Sub- 
energy solution scattering in the same range of powers has been initiated in [13]. For the 
L^-critical problem, scattering is a very hard problem, but it seems that the problem has 
now been resolved in a series of new papers by B. Dodson. For mass sub-critical solutions, 
scattering even in the energy space is a very hard problem, and is probably false. Neverthe- 
less, two particle Morawetz estimates have been used for the problem of the existence (but 
not uniqueness) of the wave operator for mass subcritical problems, |29j . We have already 
mentioned their implementation to the hard problem of energy critical solutions in [2] , [26] , 
and [11]. Recent preprints have used these inequalities for the mass critical problem, |15j . 
and the energy super-critical problem, [30J. For a frequency localized one particle Morawetz 
inequality and its application to the scattering problem for the mass-critical equation with 
radial data see [l5] . 

1.2. The Gross-Pitaevskii Hierarchy. The q-GY' hierarchy appears naturally in the 
derivation, based on the BBGKY hierarchy of density matrices, of the nonlinear Schrodinger 
equation as the dynamical mean field limit of the manybody quantum dynamics of Bose 
gases with (1 + |)-particle interactions. Recently, this research area has experienced some 
remarkable progress, see [T71 [181 UHl [SD [32l [39] and the references therein, and also [H [T6| 
[201 [211 [221 [281 HI]. 

Of particular importance for our work is the method developed in works of Erdos, Schlein, 
and Yau, [171 [HI [l9], which consists of the following main steps: 
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(1) One determines the BBGKY hierarchy of marginal density matrices for particle 
number A^, and derives the Gross-Pitaevskii (GP) hierarchy in the limit N — t- oo, 
in a scaling where the particle interaction potential tends to a delta distribution; 
see also [HTlliO]. 

(2) One establishes uniqueness of solutions to the GP hierarchy. This is the most 
difficult part of this analysis, and is obtained in [n\ [TBI [T9] by use of Feynman 
graph expansion methods inspired by quantum field theory. It is subsequently 
verified that for factorized initial data 

k 

i=i 

the solution of the GP hierarchy remains factorized for all i € I C R, 

k 

l^'"Ht,Xk;x!k) = Yl4>{t,Xj)(l){t,x'-), 
if (j){t) e H^{R'^) solves the defocusing NLS, 

for t G I C M, and (j)(0) = (po £ H'^(R'^) with q £ {2,4}. In other words, the 
solutions of the GP hierarchy are governed by a cubic NLS for systems with 2-body 
interactions {q = 2), [T71 [IHl [IHl [H], and quintic NLS for systems with 3-body 
interactions (g = 4), [3]. 

More recently, Klainerman and Machedon have introduced an alternative method in [32] 
to prove the uniqueness of solutions to the GP hierarchy for q = 2 and d = 3, in spaces 
defined via spacetime bounds on the density matrices in the GP hierarchy; these spaces 
are different from the ones used in \17\ \TE[ I19j . The uniqueness result in [32] relies on 
the assumption of an a priori spacetime bound on the density matrices. For dimensions 
d < 2, this assumption can be directly verified by exploiting energy conservation in the 
BBGKY hierarchy in the limit N — )• oo, combined with a Sobolev type inequality for density 
matrices. This was recognized in the paper [31] by Kirkpatrick, Schlein and Staffilani where 
the authors prove uniqueness for the cubic case in d = 2, and establish the assumption made 
in [32]. In [3], the corresponding problem in d = 1,2 is solved for the quintic case. Both 
|31j and [1] involve a step where a certain spacetime norm is controlled by using a Sobolev 
inequality for density matrices, and an a priori energy bound, without exploiting the decay 
in time. This approach is applicable for d < 2, but not for dimension d = 3. In d = 3, 
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a method is necessary that is truly based on spacetime norms; and such an approach has 
recently been developed in |8]. 

In [5], Chen and Pavlovic initiated the study of the well-posedness of the Cauchy problem 
for GP hierarchies with both focusing and defocusing interactions, starting directly on 
the level of GP hierarchies, and independent of the derivation from A^-body Schrodinger 
dynamics. Accordingly, the corresponding GP hierarchies are referred to as cubic {q = 2), 
quintic {q = 4), focusing, or defocusing GP hierarchies, depending on the type of the 
NLS governing the solutions obtained from factorized initial conditions. We note that for 
instance, it is currently not known how to rigorously derive a GP hierarchy from the N oo 
limit of a BBGKY hierarchy with L^-supercritical, attractive interactions. 

For the analysis of the Cauchy problem for g-GP hierarchies in [5], the following Banach 
space of sequences of fc-particle marginal density matrices is introduced 

e = {T = (7('=)(xi,...,Xfc;x'i,...,4))fceN|Tr7('=^ < oo } . (1.10) 

Given ^ > 0, the space 

K = {riliriiw- < oo} (1.11) 

is endowed with the norm 

fceN 

where 
with 

is the norm considered in [32]. If P € Ti'^, then an upper bound on the typical H°'- 
energy per particle. The parameter a determines the regularity of the solution (e.g. for 
cubic GP when d = 3, a € [l,oo)). In [5], the existence and uniqueness of solutions for 
energy subcritical focusing and defocusing cubic and quintic GP hierarchies is proven in a 
subspace of H'^. The proof involves a Picard fixed point argument, and holds for various 
dimensions d, without any requirement on factorization. The solutions derived in [5] satisfy 
a space-time bound corresponding to the one used (as an a priori assumption) in the work 
of Klainerman and Machedon, [32J. The parameter ^ > is determined by the initial 
condition, and it sets the energy scale of the given Cauchy problem. 
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The work [5] motivated further study of the Cauchy problem for q-GP hierarchy: 

• In ^9j, we identified a conserved energy functional Ei{T{t)) = i?i(r(0)) describing 
the average energy per particle, and we proved virial identities for solutions of GP 
hierarchies. In particular, we use these ingredients to prove that for L^-critical and 
supercritical focusing GP hierarchies, blowup occurs whenever £'i(ro) < 0, and 
when the variance is finite. We note that prior to [9j, no exact conserved energy 
functional on the level of the GP hierarchy was identified in any of the previous 
works, including [31] and [171 [18]. 

• In [B], new higher order energy functionals were introduced, and proven to be con- 
served quantitites for energy subcritical defocusing, and subcritical (de)focusing 
GP hierarchies, in spaces similar to those used by Erdos, Schlein and Yau in [171ll8j. 
By use of this tool, a priori bounds for positive semidefinite solutions are proven 
in the same spaces. Moreover, global well-posedness was obtained for positive semi- 
definite solutions in the spaces studied in the works of Klainerman and Machedon, 
[32], and in [5]. 

• A new proof of existence of solutions to cubic and quintic GP hierarchies for focusing 
and defocusing interactions was obtained in [7j. It does not require the a priori 
bound on the spacetime norms, which was introduced in the work of Klainerman 
and Machedon, [32], and used in the earlier work of Chen and Pavlovic [5]. 

• Very recently, a new derivation of the defocusing cubic GP hierarchy has been ob- 
tained in [8] for dimensions d = 2,3, which crucially involves generalizations of the 
tools developed in the above mentioned paper [7j. In particular, it is established 
in [8] that the GP hierarchy obtained from the N oo limit of the correspond- 
ing BBGKY hierarchy, is contained in the space introduced by Klainerman and 
Machedon in [32j based on their spacetime norms. While these results do not as- 
sume factorization of the solutions, consideration of the special case of factorized 
solutions yields a new derivation of the cubic, defocusing NLS in d = 2, 3. 

In this paper, we extend the study of Morawetz type identities for the GP hierarchy, 
inspired by our proof of a virial identity for GP hierarchies in [U] . More precisely, we derive 
an interaction Morawetz identity, and as a consequence, a regular one-particle Morawetz 
identity for solutions of the GP hierarchy. These calculations are carried out in Section [3] 
of the paper. 
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2. The nonlinear Schrodinger equation. 
We start with the equation 

iut + Au = X\u\P-^u (2.1) 

with p > 1 and A G M. We use Einstein's summation convention throughout this sec- 
tion. According to this convention, when an index variable appears twice in a single term, 
once in an upper (superscript) and once in a lower (subscript) position, it implies that we 
are summing over all of its possible values. We will also write Vj-n for For a func- 
tion a(x, y) defined on R*^ xM" we define '^xj cl{x, y) = ^'g^'^'* and similarly for Vx^k cl{x, y). 



We define the mass density p and the momentum vector p, by the relations 

P = hp, Pk = Q{uS/ku). 

It is well known, [3], that smooth solutions to the semilinear Schrodinger equation satisfy 
mass and momentum conservation. The local conservation of mass reads 

dtp + 2divp = dtp + 2Vjp' = (2.2) 

and the local momentum conservation is 

dtpi + v\5i{-\Ap + X?-^^\ur^) +ai) =Q (2.3) 

where the symmetric tensor ajk is given by 

Ujk = juVku). 

Notice that the term A^^|ti|^''''^ is the only nonlinear term that appears in the expression. 
One can express the local conservation laws purely in terms of the mass density p and the 
momentum p if we write 

A^I.r^ = 2^A^p^ 
p+l p+l 

and 

ajk = 2^{S/juVku) = -{2pjPk + ^'^jP^kP), 
P ^ 

but we won't use this formulation in this paper. We are ready to state the main theorem 
of this section: 

Theorem 2.1. UHl [131 ESI E] Consider u G Ct(M;C^(R")) a smooth and compactly 
supported solution to (12. ip with u{x,0) = u{x) € C^(M"). Then for n>2 we have that 

^,1^ n-3,, ,o,,,9 , ,x^P— 1 f f \u(y,t)\'^\u(x,t)\'^^^ , , , 

C\\D — ^{\u\^)\\}2r2 + {n-l)X- / / ' , ' ^ dxdydt 



2,2^2 -r V - — --r 



\x - y\ 
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< ||«o|li2 sup \My{t)\, 

tGM 

where 

My{t) = [ ^"^ • 'i{u{x)Vu{x))dx, 

Jk" f ~ y\ 

D°' is defined on the Fourier side as D°^f{£^) = for any a G M and C is a positive 

constant that depends only on n, j42j . For n = 1 the estimate is 

\\dx{\u\^)\\l2L2 +A^^-^||u||^++3 +3 < ^i|uo||i2 sup\\d^u\\L2. 

Remarks on Theorem 12.11 

1. By the Cauchy-Schwarz inequality it follows that for any n >2, 

sup\My{t)\ < ||no|||2 sup ||Vn(t)|li2. 
A variant of Hardy's inequality gives 

sup!M^(t)| < ||no|||2 sup 1 , 

For details, see |24j . 

2. Concerning our main theorem, we note that both the integrated functions in the second 
term on the left hand side of the inequalities are positive. Thus when A > 0, which 
corresponds to the defocusing case, and for data say, we obtain for n > 2: 

P"'^(kl')llL?Li < holli snp\\Vumh<Miuof^E{uo)K 

tm 

and for n = 1 

\\dxi\u\^)\\l2L2 < \\uo\\l,snp\\dMt)\\l2<Miuo)lE{uo)l 

' ^ t€R 

These are easy consequences of the conservation laws of mass (I1.6P and energy (jl.Sp . They 
provide the global a priori estimates that are used in quantum scattering in the energy 
space, [23]. 

3. Analogous estimates hold for the case of the Hartree equation iut + Au = A(|x[^'''* \u\'^)u 
when < 7 < n, n > 2. For the details, see [29]. We should point out that for < 7 < 1 
scattering fails for the Hartree equation, [27], and thus the estimates given in }29] for n >2 
cover all the interesting cases. We also expect the Hartree quadratic Morawetz estimates 
that were established in [29] to extend to the case of the Hartree hierarchy, for which 
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factorized solutions are determined by the Hartree equation. A detailed analysis of this 
problem for GP hierarchies is presented in section [3l below. 



4. Take A > 0. The expression 

\\D-^{\u\^)\\^Ll^ 

for n = 3, provides an estimate for the L^L^ norm of the solution. For n = 2 by Sobolev 
embedding one has that 

For n > 4 the power of the D operator is negative but some harmonic analysis and inter- 
polation with the trivial inequality 

ll-D5n||ioci2 < ||n|| .1 
provides an estimate in a Strichartz norm. For the details see |44] . 

5. In the defocusing case all the estimates above give a priori information for the Hi- 
critical Strichartz norm. We remind the reader that the iJ'^-critical Strichartz norm is 
||u||^9^r where the pair {q,r) satisfies | + 7 = § ~ s. In principle the correlation of k 
particles will provide a priori information for the H^^k critical Strichartz norm. In Id an 
estimate that provides a bound on the H» critical Strichartz norm has been given in |12j . 

6. To make our presentation easier we considered smooth solutions of the NLS equation. 
To obtain the estimates in Theorem 12.11 for arbitrary functions we have to regularize 
the solutions and then take a limit. The process is described in [24] . 

7. A more general bilinear estimate can be proved if one correlates two different solutions 
(thus considering different density functions pi and ^2)- Unfortunately, one can obtain 
useful estimates only for n > 3. The proof is based on the fact that — A^[x| is a positive 
distribution only for n > 3. For details the reader can check [13]. Our proof shows that the 
diagonal case when pi = p2 = |up provides useful monotonicity formulas in all dimensions. 

Proof. We define the Morawetz action centered at zero by 

Mo(t) = / Va(x) dx, (2.4) 

where the weight function a{x) : ^ M is for the moment arbitrary. The minimal 
requirements on a{x) call for the matrix of the second partial derivatives djdka{x) to be 
positive definite. Throughout our paper we will take a{x) = \x\, but many estimates can 
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be given with different weight functions, see for example [llj and If we differentiate 
the Morawetz action with respect to time we obtain: 



dtMo{t) = / Va(x) • dtp{x) dx = V ja{x)dtp' (x) dx 

= / {VjV^a{x))5l{-]-Ap + X^^^\u\P+^)dx + 2 [ {VjV^aix))^V^uVku)dx, 

where we use equation (j2.3p . We rewrite and name the equation as follows 

dtMo{t)= [ Aa{x){-lAp + X?^\u\P+^)dx + 2 [ {VjV^a{x))^V^uVkU)dx. (2.5) 

Notice that for a(x) = \x\ the matrix VjVfca(x) is positive definite and the same is true if 
we translate the weight function by any point y € M" and consider Va;jV^''^a(a; — y) for 
example. That is for any vector function on M", with values on M or C we have 

that 

{VjV''a{x))v^{x)vkix)dx > 0. 



To see this, observe that for n > 2 we have Vj-a = and VjVka = j^j (^Skj — 



Summing over j = k we obtain Aa{x) = ^^j^- Then 

VjV'a{xyix)v,ix) = ^{6k, - ^)v^ix)v,ix) = ^(\vix)\' - (^^)') > 

\x\ kc r \x\ \ \x\ ' / 



by the Cauchy-Schwarz inequality. Notice that it does not matter if the vector function is 
real or complex valued for this inequality to be true. In dimension one (12. Sp simplifies to 



dtM^{t)= [ a,,{x){-lAp + X^-^\u\P+^ + 2\u,\^)dx. (2.6) 
Jr 2 p+1 

In this case for a(x) = we have that axx{x) = 25{x). Since the identity (12. 5p does not 

change if we translate the weight function by y G M" we can define the Morawetz action 

with center at y G M" by 



My{t) = / Va{x - y) ■ p{x) dx. 
We can then obtain like before 

dtMy{t)= [ Axa{x-y){-lAp + X^\u\P+')dx 

+2 / (V^jV^''=a(x-y))5?(V^'%V^,fc'u)drc. 
We now define the two-particle Morawetz action 

M(t) = I \uiy)\^Myit) dx 
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and differentiate witfi respect to time. Using tlie identity above and tlie local conservation 
of mass law we obtain four terms 

dtM{t) = [ \u{y)fdtMy{t) dx+ [ dtp{y)My{t) dx 

= [ \u{y)\^AMx-y){-lAp + X^\uf+')dxdy 

+2 / |t/(2/)|2(V^jV"'*^a(x-y))$R(V^'%V^,fc?i)dx(iy 

-2 / Vy'^pj{y)V^,ka{x - y)p''{x)dxdy 

= 1 + 11 + 111 + 2 [ pj (y) V^'-J V^,ita(x - y)p^{x)dxdy 
by integration by parts with respect to the y— variable. Since 

V^'^ V^,ika(x -y) = -V^'^V^,fea(a; - y) 

we obtain that 

dtM{t) = 1 + 11 + 111-2 f V^'^V,,,fca(x - y)pj{y)p^{x)dxdy (2.7) 



= 1 + 11 + 111 + IV 



where 



1=1 \u{y)fA^a{x-y){-l-Ap)dxdy, 
11= f \u{y)fA,a{x-y){\^\u\P+')dxdy, 
111 = 2 [ |u(y)|2(V^,,V^''=a(x - y))5?(V^'%V^,feu)dxdy, 
IV = -2 I V^'^ V^,jka(a; - y)pj{y)p\x)dxdy. 

^K^xRj 

Claim: 77/ + IV > 0. Assume the claim. Since Aj.a{x — y) = ^E^\ we have that 

d^M{t) > ^ / Ap)dxdy + (n - 1)A^ / Myl^|^,(^)|P+id;^dy. 

2 j]R!}xMy F ~ y| J* + 1 JMrfxR? F ~ y\ 

But recall that on one hand we have that — A = and on the other that the distributional 
Fourier transform of |^ for any n > 2 is where c is a positive constant depending 

only on n. Thus we can define 

D-(--')f{x)=C [ ji^dy 

^R" F - y\ 



i2|i2 
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and express the first term as 

2 JMyxIRJ \X — y\ 2 

by the usual properties of the Fourier transform for positive and real functions. Integrating 
from to t we obtain the theorem in the case that n > 2. 

Proof of the claim: Notice that 
III + IV = 2 [ V^jV^'^a(x - y)(|u(y)|2sft(V^'%(x)V^,fetx(x)) - {y)pk{x)) dxdy 

= 2 / V,jV^''^a(x-2/)(44^(^(^)(V"'M^)M^)(V.,itu(x))) - (y)vk{x))dxdy. 
Since 

^x,j^x,k^{x -y) = Vy,j'^y,ka{y - x) 
by exchanging the roles of x and y we obtain the same inequality and thus 

III + IV= [ V^jV^M^ - y) (44^(^(^)(V"'%(^))n(x)(V,,fcu(x))) - jP{y)pk{x) 

+^5}?(n(y)(V^'%(y)My)(Vj,,fen(2/))) - {x)pk{y)) dxdy . 
Now set zi = u{x)V x,k'U'{x) and Z2 = u{x)'V^'^u{x) and apply the identity 

to obtain 

3?(t^(x)(V^'%(x))iZ(x)(V^,fe«(x))) =3ff(«(x)V^,feu(x))K(u(x)V^'%(x)) 

+'^{u{x)Vx,ku{x))Q{u{x)W''^^u{x)) = ^Vx,kpix)V'''^pix)-pk{x)p'{x) 
and similarly 

^uiy){vyMy)My){^yMy))) = \'^yMy)^''''piy) -Pk{y)p^{y)- 

III + IV = ] [ V,jV"'*^a(x - y)44v.,ikP(^)V^'V(x)dxdy 

^jR^jxRg P{y) 

+ I S7yjVy''a{x-y)(^pk{x)pi{x)+^Pk{y)p^{y)-Pk{x)p^{y)-^^^^ 

j» ^p{x) ply) / 



Thus 



.yxRS 
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Since the matrix V x^jV'^'^a{x — y) = VyjVy'^a{x — y) is positive definite, the first two 
integrals are positive. Thus, 

/// + IV > 

[ J V'^'^a(x -y)( ^^Pk {x)p> (x) + ^^Pk {y)p' (y) -Pk {x)p^ (y) -Pk {y)p' {x)] dxdy. 

Now if we define the two point vector 



T I \ \ p{y) I \ \ p{x) , . 



we obtain that 



/// + IV > [ V^jV^''^a(a; - y)J^{x,y)Jk{x,y)dxdy > 



y ""^a; 



and we are done. 

The proof when n = 1 is easier. First, an easy computation shows that if a{x, y) = \x — y\ 
then dxxO-ix-iy) = 25{x — y). In this case from ()2.7p we obtain 

dtM{t) = I \u{y)\^26ix-y){-lpxx)dxdy + 2 [ \u{x)\^ {X^\u{x)\P+')dx 

+4 / \u{xy(^\unS^ dx - \ \ p^{x)dx. 
JR Jr 

/ \u{yy(^2b(x - y){-\-px^dxdy = [ (dx\u{x)\A dx. 

jR-uXR-^ ^ JR ^ ' 



But 



In addition a simple calculation shows that 



Thus 

4|n(x)l2lu,|2-4p2(x)= 

and the identity becomes 

dtM(t) = 2[ (dJu\'^ydx + 2 [ \u(x)\'^(X^^^\u(xW+^)dx (2.8) 
Jr^ ' Jr ^ p+l 

which finishes the proof of the theorem. □ 
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3. MORAWETZ IDENTITIES FOR THE GP HIERARCHY 

In this section, we derive one-particle Morawetz inequalities for GP hierarchies in Theo- 
rem [3?T1 and interaction Morawetz identities for GP hierarchies in Theorem 13.2^ below, for 
1 < n < 3 dimensions. For simplicity of exposition, we will only present the case of cubic 
GP hierarchies here, that is, g-GP hierarchies with q = 2. The case of quintic or higher 
degree g'-GP hierarchies > 4, g G 2N) can be treated in a completely analogous manner. 

For convenience of exposition, we assume that F = (7^^^^) solves the cubic GP hierarchy 
()1.2p in Ti^ (the corresponding local well-posedness theory is covered in [5]), with interaction 
terms defined in (jl.3p with g = 2. 

3.1. One-particle Morawetz identities for the cubic GP hierarchy. The density 
function corresponding to the one-particle marginal is defined by 

p(x) := 7(^)(x;x). (3.1) 

One can straightforwardly verify that 

dtpix) = [ dudu'e'^''-'''^^da'-^\u;u') 



- [ dudu' e^'^'^-'^'^^u'^ - {u'f)^^^\u;u') 



+ - / dudu'e^("-"')^5i,27(2)(n;n') (3.2) 



I 



= i j dudu' e^^''-'''y-'{u + u'){u-u')^^^\u;u') 

= -V^. •ydndn'e^("-">(n + n')7^^)(u;n'), (3.3) 
so that in analogy with (|2.2p . 

9tp(x) + 2V^ • P(a;) = , (3.4) 

with 

P{x) := ydndn'e*("-">^^^7W(n;n') 



(3.5) 



is the momentum operator, see also [9] . The fact that the interaction term ()3.2p equals zero 
is proven in eqs. (5.5) - (5.8) in [9]. We are here adopting conventions analogous to those 
in the previous chapter applied to the NLS. In the sequel, we will suppress the dependence 
on t from the notation, for simplicity. 



^As a result, all expressions in the steps below can easily be seen to be well-defined. By adopting the 
arguments of [S], it in fact suffices to consider solutions in T-L^. 
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We define the Morawetz action 

Ma := j dxVa{x) ■ P{x) (3.6) 
in analogy to (j'2.4|) . The time derivative is given by 

dtMa = j dxVa{x)-dtP{x). (3.7) 
Then, we obtain the following version of the regular Morawetz identity. 

Theorem 3.1. Under the conditions formulated above, the one-particle Morawetz identity 

dtMa = -\ j dx{^^^^a{x))-i^^\x;x) 

+ ^y" (A^a(x))7(^^(x,x;x,x) 

+2 3? /" dxdx'5{x - x') {d^^d^^aix))d^^d^p^^\x;x') . (3.8) 

holds for solutions of the cubic GP hierarchy. 

We defer the proof to section 13.31 
3.1.1. Factorized solutions. Substituting factorized solutions of the form 

k 

l^''Ht,Xk,x!;,) = ll<p{t,x,mt,x'^), (3.9) 
i=i 

where 

idt(j)it,x) + Aa;0(i,x) = X\<j){t,x)\'^ (j){t,x) (3.10) 
with initial data 0(0, ■) = (j)o G H^, we obtain the following result. 
dtMa = J dx{A.,AMx))m,x)\^ 

+ ^ I dx{A^a{x))\(t){t,x)\^ 



-2^jdxY^ {d,^d,,a{x)){d^,cP{t,x)) d,'^cl,{t,x') . (3.11) 



This corresponds to the one-particle Morawetz identity (12. 5p for the NLS, in the cubic case 
p = 3. 

3.2. Interaction Morawetz identities for the cubic GP hierarchy. In this section, 
we derive interaction Morawetz identities for GP hierarchies which generalize those for the 
NLS. 
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3.2.1. Morawetz action. We write 

- '.-»v,<^)(„,„;<„0. (3.12) 

For a function a : M'* x M" — >■ M, (x, y) a{x, y), we define the Morawetz action 

Ma := J dxdyVxa{x,y) ■ Px{x,y) . (3.13) 

where 

Px{x,y) := j duduUvdv' e'^''-'''^''+'^-"-'''^y (^^-^^^ (3.14) 

so that 

P{x) = j dyPx{x,y). (3.15) 

The time derivative is given by 

dtMa = J dxdyVxa{x,y) ■ dtPx{x,y) 

= i[(Ai) + (A2) + (A3) + (A4)], (3.16) 

where the four terms on the rhs are defined as follows. 
We have 

{Ai) := j dxdy{VMx,y)- j dudu' dvdv' e^^''-''>+'^'"-'''^y 

{u + u') {u^ - {u'f) 7(2) (u, V- u', v') (3.17) 

and 

(^2) := ^ Jdxdy{ V^a(x, y) ■ J du du' dv dv' e^("-"')^+i(^-^'')2/ 

{u + u'){v'^ - {v')'^)j^^\u,v;u',v'). (3.18) 

Moreover, 

{As) := ^ J dxdy{Vxa{x,y)- J dudu'dvdv'e'^''-''">''+'^''-'''^y 

{u + u') (-61^37(3) (u, V- u', v') - Si"37(3) (u, v; u' , v')) (3.19) 

and 

(A4) := ^ j dxdy{Vxa{x,y)- j dudu' dvdv' e'^''-''>+'^''-'''^y 

{u + u') (S+3^3) (^^ _ 5^3) (^^ ^. ^/)) . (3.20) 

We now discuss each of these four terms in detail. 
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3.2.2. The term (Ai). We have 
(^i) 

= ^ J dxdy V:ca{x, y) ■ J du du' dv dv' e^(«-«')^+i(''-^')2/ + u') {u^ - {u'f) 

= - j dx dyVxdixjy) ■ J dudu' dv dv' ^^'^\u,v;u' ,v') 

[{u + u') (g) (n + u')]{u - u') e^(«-«')^+i{"-^')2/ 
= -jdxdy jdudu'dvdv'^i-\u,v,u',v') 

VMx,y) ■ [{u + u') ®{u + n')](V^e^("-"')^+^(''-"')^ ) 
= JdxdyJdudu'dvdv-,^^Hu,v;u',v') (3.21) 

where u := u + u' . This equals 

= j dxdx' dy 5{x — x') j dudu'dvdv' y^^{dxidx^a{x,y)) 



1,3 



i{ux-u'x')\i{v-v')y ^(2) ^. ^/^ ^/-) 



- / dxdx' dyS{x - x') ^{dxidx^x^y)) {d^i - dx^){dx^ - dx'.)'j^^\x,y;x' ,y) 

- J dxdx' dy 6{x - x') Y^{dx,dxja{x,y)) {d^.d^j + dx'.dx'.)'y^'^\x,y;x' ,y) 
+ / dxdx' dy S{x - x') '^{dx,dxja{x,y)) {d^'.dxj + d^^^^dx'^) 'y^^\x,y;x' ,y) 

- / dxdx' dy 6{x - x') ^{dx^dx^x.y)) 

(dxidxj + dxidx'. + dx'.dxj + dx'.dx'. ) 7^^^ (x, y, x', y) 
+ 2 / dxdx' dy 5{x - x') ^{dxidx^x^y)) {dx'.dx^ + dxidx'.)-i'''^\x,y;x' ,y) 

dxdy ^{dxidxja(x, y)) dxidx^ 7^^^ ix, y; x, y) 

i,j 

+ mj dxdx'dy5{x-x') Y,{dx,dx,aix,y))dx,dx'.7^^\x,y;x',y) (3.22) 
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where we used ^^"^^x , y; x' , y) = j^'^\x',y;x,y), and applied a coordinate change x -H- x' in 
one of the two integrals contributing to the last line. This equals 



dxdx' dy (A^a(j;, y)) A^.7(^) (x, y; x, y) 
+ m j dxdx'dy5{x-x') ^(a,,a,.^,a(x,2/))9,,9,/7(2)(x,y;x',y). (3.23) 

This corresponds to the first and last term on the rhs of (3.36) in [43J . 
3.2.3. The term {A2). We have 

(^2) := ^ Jdxdy{ V^a(x, y)) ■ j du du' dv dv' e^(«~«')^+i{^-^')!/ 
[(n + u') (8> (f + v')]{v — v') 7^^^ {u, v; u' ,v') 
= - I dxdy^{ d^^a{x, y)) ■ I du du' dv dv' (ay^e'("-"')^+*(''-^')2') 

{u + u')j{v + f')^7*-^^(u, u ,v') 
dx dy ( d^^dy,a{x, y) ) (3.24) 

du du' dv dv' f,i{u-u')x+i{v-v')y ^ ^ ^/^^ ^(2) ^. ^/^ ^/^ _ 



Here we note that if j^'^\x,y;x',y') = j^^\x;x')j^^\y;y') has product form, then the 
integral on the last line corresponds to 4P(x)P(?/), the product of momentum densities 
defined in (13.511. 



3.2.4. The term {A3). We have 

5+37(3) (x,y;x',y') 

dz dz' 6{x — z) 5{x — z') 

j dudvdqdu'dv'dq'e''^'''-'+^y+'''-''''''-'''y'-'^'''^ 7(2) (u, v, q; u',v', q') 

^ dudvdqdu'dv'dq'e'^^''+''-''">''-''''''+^y-^'y'^ 7(2) (n, v, q; u',v', q') . (3.25) 
Therefore, 

B^^){u,v-u',v') 

dxdx'dydy' ^-iux-ivy+iu'x'+iv'y' (5+^^(3))(^^ y. ^> ^ y/) 

dqdq t'^'^-' {u — q + g', f , q; u ^v' , q) . (3.26) 
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Likewise, one obtains 

Bf^^){u,v;u',v') = Jdqdq'^(^\u,v,q;u' + q-q',v',q'). (3.27) 
Now, in order to consider {As) we first look at 

^ J dudu'dvdv' e^(«-"')a^+i{''-^')2/ {u + u') {B+~f^) {u, v; u' , v') - B^^) {u, v; u' , v')) 
= \ j dudu'dvdv' dqdq' e^(«-«')^+i(^-''')2/ {u + u') 7^=^) {u-q + q', v, q; u' , v', q') 

J dudu'dvdv' dqdq' e*("-"')^+i(^-'^')2/ {u + u') 7(2) {u, v,q;u' + q- q', v', q') .(3.28) 

In the last term, we apply the change of variables u ^ u — q + q' and u' ^ u' — q + q', so 
that the difference u — u' remains unchanged. We obtain that the above equals 

^ j dudu'dvdv' dqdq' e^(«-«')a^+i{''-^')2/ {u + u') 7^^) {u-q + q', v, q; u',v', q') 

j dvdv'dqdq' e^("-"')^+i(^-''')2/ (^u + u'-2q + 2q') 7^^) {u-q + q', v, q; u' , v', q') 
= \ j dudu'dvdv' dqdq' e*('^-'^')a;+i(^-^'')2/ ( + _ ^u + u' -2q + 2q') ) 

7^^) {u-q + q', v, q; u', v', q') 

= - f dvdv'dqdq' e'(n-u')x+i{v-v')y _ g') ^(3) {u - q + q' ,v,q;u' ,v' ,q') . (3.29) 



The contribution of this term to {A3) is given by 



j dxdyV^a{x,y)- j dudu'dvdv'dqdq' 6^^''-'''^''+'^''-'''^^ 

2{q - q') 7(3) {u-q + q' , v, q; u' , v' , q') . (3.30) 



Next, we express everything in position space. 
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We have that the last hne equals 
^ J dxdy V^a(x, y) ■ J dXdYdZdX'dY'dZ' j dudu' dvdv' dqdq' (^^u-u')x+{v-v')y _ g') 

^i{-{u-q+q')X-vY-qZ+u'X'+v'Y'+q'Z')^(3) (^j^ ^ y, Z ; X' ,Y' , Z') 

= ^ J dxdy j dXdYdZdX'dY'dZ'j^^\X,Y,Z;X',Y',Z') j dudu' dvdv' dqdq' 

^iuix-X)+iv{y-Y)~iu'(x-X')-iv'{y-Y') 2V.^a{x,y) ■ {q - q') e+i<J{X-Z)-q'{X-Z') 

= j dxdy j dXdYdZdX'dY'dZ' -f^^\X,Y,Z;X',Y',Z') j dqdq' 

6{x - X)d{x - X') 6{y - Y)6{y - Y') 2Vxa{X, Y) ■ Vx e+^9(^-^)-^9'(^-^') 

= j dXdYdZdZ'-f^^\x,Y,Z;X,Y,Z') 

2Vxa{X,Y) -VxSiX - Z)6{Z - Z') (3.31) 

= -A y dXdYdZ 7(2) (X, Y, Z; X, Y, Z) 2Vxa(X, Y) ■ Vx S{X - Z) 

= A y dXdYdZ 6iX - Z) (2Axa(X, Y) + 2(Vxa(X, Y)) ■ Vx ) 

7(3)(X,y,Z;X,y,Z) (3.32) 

where we have written 5{X - Z)5{X - Z') = 5{X - Z)6{Z - Z') to get dOB . 
Now we rename the variables {X, Y, Z) —?■ (x, y, z), and note that 

J dxdy {Vr,a{x, y)) ■ V^^i^^'' (x, y, x; x, y, x) 
= j dxdydz6{x - z) ( (Vx.a(x, y)) ■ + (Vza^z, y)) ■ )7*^^) {x, y, z; x, y, z) 

= j dxdydz6{x - z){{Vxa{x,y)) ■Vx'y^^Hx,y,z;x,y,z) 

+ C^zaiz, y)) ■ V27(^) {z, y, x; z, y, x) ) 
= J dxdydzS{x - z){2{Vr,a{x,y)) ■VxJ^^\x,y,z;x,y,z)) (3.33) 

where we used the symmetry 'y^^\x,y, z; x,y, z) = j^^\z,y,x; z,y,x), and renamed the 
variables in the last term. Clearly, the left hand side equals 

- J dxdy{Axa{x,y))'y^^\x,y,x;x,y,x) (3.34) 

from integrating by parts. 
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Therefore, combining I^M), ([3:33]) and (l33i]) 

{A3) = X dxdydz6{x - z){2Ar,a{x,y) - A:^a{x,y))-y^'^'>{x,y,z;x,y,z) 

This corresponds to the second term on the rhs of (3.36) in |43j . 

3.2.5. The term (A4). We have 

(^4) = 7/ dxdyV^a{x,y)- J dn' di- dt;' e^("-"')^+*(''-^')s' 

{u + u') [Bl^3)(u,v;u',v') - B^^3){u,v;u',v')] 

=■■ ^ j dxdy\7^a{x,y) ■ j dudu' dv dv' dqdq' e'^''-'''>+'^''~'''^y (3.36) 
{u + u') [^^^^UjV - q + q' ,q;u' ,v' ,q') - j^^^u,v,q;u' ,v' + q - q' ,q')] 

= ■■ \ j dxdyV^a{x,y) ■ J du du' dv dv' dq dq' e'^''-''">'=+'^''-'"">y (3.37) 

{u + u') [7(3) (n, V + q', q; u',v' + q, q') - 7(2) (n, v + q' , q; u',v' + q, q')] 
= (3.38) 

where to pass to (I3.37p . we used the coordinate change v ^ v + q, v' ^ v' + q for the 
expression involving the first term in the square bracket in (j3.36p . and v v+q' , v' — ?> v'+q' 
for the second term. Both coordinate changes leave the difference v — v' invariant. 

3.2.6. Completing the proof. Summarizing, we obtain from p.l6p that the following result 
holds. 

Theorem 3.2. Under the conditions formulated above, the interaction Morawetz identity 
dtMa = -\j dxdy{Axa{x,y))Ax'y^'^\x,y;x,y) 

+ ^ I dxdy{Ar,a{x,y))-f^^\x,y,x;x,y,x) 



+ 23f?y" dxdx'dy5{x - x') ^ {dx.dx^a{x,y))d;red^p^^\x,y;x' ,y) 
+ 2 j dxdyY, [d^.dyA^^v)) j dn dn' d^; d^;' e^("-"')^'+*(''-^')^ 



' u + u'\ f V + v' 



2 V 2 



^^"^Xu^v^u^v') . (3.39) 



holds for solutions of cubic GP hierarchies. 
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We compare this result with ()2.7p obtained above for the NLS. 

For factorized solutions of the cubic GP hierarchy of the form ()3.9p . we obtain the 
following: 

dtMa = j dxdy\(t>{y)\'^{AMx,y))AxP{t,x) 

+ dx{AMx,y))m,y)\^m,x)\^ 

+ 25? / dxdy\cl){y)\^ ^ ( d^^d^,a{x, y) ) x)) {8^^ <p{t, x)) 

+ 2 j dxdyY, {dx, dy^a{x, y) ) {P^t, x))j {P^{t,y))i , (3.40) 



where 



P^x) := ^i<j,{x)V.<f>{x)) (3.41) 

is the momentum density corresponding to <f>{t,x). This corresponds to (j2.7p for the cubic 
NLS where p = 3 (so that = In particular, we note that for a{x,y) = a{x — y), 

we evidently have dx^dy^a{x,y) = —dx^dx^aix^y), which agrees with ([2 



3.3. Proof of the one-particle Morawetz identities. In this section, we prove the stan- 
dard (single-particle) Morawetz identities in Theorem 13. H as a corollary of the interaction 
Morawetz identities derived above. 

The corresponding explicit expression for dtMa can be easily obtained from the interac- 
tion Morawetz identities (I3.39p . by choosing 

a{x,y) = a{x) (3.42) 

independent of y, and 

^(i)(x;x') = / dyj^^\x,y;x',y) 



-f'^'^\x,z;x',z') = / dyj^^\x,y,z;x',y,z'), (3.43) 



which follows from the admissibility of the density matrices, see (|1.4p . 
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Accordingly, ()3.39p reduces to 

X X J X ^ X J 

+2^ I dxdx'5{x - x') {d^^d^,aix))d^,d^>^-f^^\x;x') . (3.44) 

We note that the term involving the momentum densities on the last line of ()3.40p is not 
present here (since dy^a{x) = 0). 
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